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ABSTRACT. We review the use of peculiar velocities of galaxies as a probe of 
cosmological models. We put particular emphasis on comparison of the pecu- 
liar velocity and density fields, focussing on the discrepancies between various 
recent analyses. We discuss the limitations of the commonly used linear bias 
model, which may lie at the heart of some of the current controversies in the 
field. 



1 Introduction 

We believe for a large number of rea- 
sons (including those to be outlined 
in this review) that we live in a uni- 
verse which is dominated by dark 
matter. We see galaxies by the light 
of the stars that shine within them; 
we can similarly infer the presence of 
large quantities of hot gas through its 
X-ray emission, and cold gas through 
molecular and atomic emission in the 
millimeter and radio portions of the 
spectrum. But the baryons that give 
rise to these emissions comprise at 
most 10% or so of the mass density 
of the universe. 

The study of large-scale structure 
asks how this dark matter is dis- 
tributed in space. We cannot observe 
it directly, but we can infer many of 
its properties due to its gravitational 
influence on the matter around it. 
We observe that the distribution of 
galaxies is clustered, and make the 
assumption (which is at least qualita- 
tively justified, as we will see below) 
that the dark matter distribution is 
clustered in a similar way. A local 
overdensity in mass will gravitation- 
ally attract the adjacent matter, dis- 
turbing the pure Hubble flow; simi- 
larly, a local underdensity will grav- 
itationally repel matter. These addi- 
tional motions, or peculiar velocities, 
are observable. Therefore, we can use 
peculiar velocities to get an observa- 
tional handle on the clustered com- 
ponent of the dark matter. 
We make these notions concrete as 
follows. 

Pure Hubble flow states that redshift 
is proportional to distance, 



and indeed, if we measure distances 
in units of km s _1 , and never make 
reference to standard yardsticks cal- 
ibrated in physical units (e.g., cm or 
Mpc), the value of Ho is identically 
unity. Peculiar velocities are defined 
relative to the comoving expanding 
frame implied by equation (hi): 



r + f ■ [v(r) - v(0) 



(2) 



where v(r) is the peculiar velocity 
field, and redshifts are measured rel- 
ative to the barycenter of the Local 
Group (as opposed to, e.g., the CMB 
rest frame). 

In linear gravitational instability 
theory one can relate the density 
fluctuations to the peculiar veloci- 
ties: 



An J | r — r | J 



where Q is the Cosmological Density 
Parameter, and the density fluctua- 
tion field is defined as: 



5(r) 



P(r) - (p) 
(p) ' 



(4) 



The integral on the right-hand side is 
just proportional to the gravity vec- 
tor. As an aside, note that it is man- 
ifestly not true that peculiar veloc- 
ity is proportional to gravity in the 
highly nonlinear regime: the gravity 
and velocity vectors of the Earth's 
orbit are at right angles to one an- 
other! 

We will also find it useful to express 
equation (^) in differential form; by 
taking the divergence of both sides of 



cz = Hgr, 



(1) 
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the equation, one finds: 

V ■ v(r) = -n os 8(r). (5) 

Note that peculiar velocities, if mea- 
sured in the linear regime, are a mea- 
sure of the response of galaxies to 
the gravitational influence of all mat- 
ter, and especially the dominant dark 
matter. Therefore, we can measure 
the large-scale distribution of dark 
matter directly from observations of 
the peculiar velocity field. 
Moreover, if we have an indepen- 
dent measurement of the large-scale 
distribution of galaxies from redshift 
surveys, and a model that relates the 
galaxy and dark matter distribution, 
we can use equation (jsj) or (jHJ) to 

• Test the gravitational instability 
paradigm from which these equa- 
tions were derived; 

• Test the assumed model relating 
the distribution of the dark and 
luminous matter; 

• Measure the quantity Q directly. 

2 The Measurement of Pecu- 
liar Velocities 

Peculiar velocities manifest them- 
selves through their modification 
to the Hubble Law, equation (Q). 
Therefore, with independent mea- 
surements of redshifts and distances 
of a galaxy, we can measure the ra- 
dial component of its peculiar veloc- 
ities. Redshifts are relatively easy to 
measure; a high-quality galaxy spec- 
trum will show sharp stellar absorp- 
tion lines, and often strong emis- 
sion lines that arise from photo- 
ionized gas in Hn regions or an ac- 
tive nucleus. Distances are quite a bit 
more difficult; one needs to identify 
a standard candle (i.e., a distance- 
independent way to determine the 
absolute magnitude of the galaxy, or 
some well-defined component of it), 
or a standard yardstick (i.e., a char- 
acteristic physical size of the galaxy). 
Comparison with the apparent mag- 
nitude or angular size then yields a 
distance. 

Reviews of distance indicators used 
in peculiar velocity work can be 
found in Jacoby et al. (1992), 
Strauss & Willick (1995, hereafter 
SW), and Willick (1998). The Tully- 
Fisher (1977: TF) relation states 



that the luminosity and (inclination- 
corrected) rotation velocity of a spi- 
ral galaxy are related by a power law; 
expressing this in terms of absolute 
magnitudes, we find: 

M = AlogAV + B. (6) 

Thus, observations of apparent mag- 
nitudes m and rotation velocities A V 
(via the Ha line at 6563A with a 
long-slit spectrograph placed along 
the major axis of the galaxy, or the 
width of the 21 cm line in the radio), 
yield a distance. The TF relation 
is not perfect, and shows some ap- 
preciable scatter; distances are mea- 
sured typically to an accuracy of 15- 
20%. Moreover, the constants A and 
B depend on the details of sam- 
ple selection, the band in which the 
galaxy photometry is done, and so 
on, and these quantities need to be 
calibrated carefully for any given TF 
sample. Since we measure the pecu- 
liar velocity as the difference of two 
rather large quantities (equation ^), 
the signal-to-noise ratio per galaxy 
of the peculiar velocity field is typ- 
ically low, and is a decreasing func- 
tion of distance. This means that a 
great deal of care needs to be taken in 
doing TF work; when data are noisy, 
there are a number of pernicious sys- 
tematic effects that can creep into 
the derived velocity field if one is 
not careful (cf., the discussion in SW; 
Teerikorpi 1997). 

There are a number of other distance 
indicators that are used in peculiar 
velocity work: 

• Elliptical galaxies are observed 
to fall along the so-called funda- 
mental plane (Djorgovski & Davis 
1987; Dressier et al. 1987), which 
relates the luminosity, the surface 
brightness, and the velocity dis- 
persion or color of galaxies. This 
yields distances to an accuracy of 
roughly 20% in distance, compa- 
rable to the TF relation. 

• Distances can also be measured 
to ellipticals using the method 
of surface brightness fluctuations 
(Tonry & Schneider 1988). Ellip- 
tical galaxies typically have in- 
trinsically smooth surface bright- 
ness profiles, unmarred by struc- 
tural components such as spiral 
arms or dust lanes. However, their 
light is generated by a finite num- 
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ber of stars, and the closer the 
galaxy is to us, the fewer of the 
red giant stars that dominate the 
optical light will fall within each 
seeing element of an image. This 
number of stars will be subject 
to Poisson statistics, and there- 
fore the closer the galaxy is to 
us, the rougher the image will 
appear. This roughness is there- 
fore a measurement of the dis- 
tance of the galaxy. This is a 
method that holds a tremendous 
amount of promise for peculiar 
velocity work, as it yields dis- 
tances to accuracies of 5% for 
at least some galaxies. Unfortu- 
nately, the method is seeing lim- 
ited, and requires the resolution 
of the Hubble Space Telescope to 
reach much beyond 3000 km s _1 . 
Although there now exist sub- 
stantial data for nearby galaxies 
using this technique (e.g., Tonry 
et al. 1997), it has not yet been 
exploited for studies of the pecu- 
liar velocity field in any detail. 

• The most luminous galaxies in 
clusters have long been recognized 
to be distinct from the "ordinary" 
ellipticals that make up the bulk 
of the galaxy population in ellip- 
ticals. Gunn & Oke (1975) and 
Postman & Lauer (1995) show 
that the luminosity of Brightest 
Cluster Galaxies within a met- 
ric aperture is directly related to 
the logarithmic slope of their sur- 
face brightness profile; using this 
yields a distance to an accuracy 
comparable to that the TF re- 
lation, 15-20%. These have been 
used to measure the bulk compo- 
nent of the velocity field on very 
large scales (Lauer & Postman 
1994); more on this below. 

• Finally, it has long been rec- 
ognized that Type la super- 
novae are standard candles (e.g., 
Phillips 1993). When one includes 
information about the shape of 
the light curve, the peak lumi- 
nosity is predicted well enough 
to yield distances accurate to 5% 
(Hamuy et al. 1996; Riess, Press 
& Kirshner 1996). Supernovae are 
of course rare events, and there- 
fore one cannot a priori pick one's 
sample of galaxies via this tech- 
nique; nevertheless, there have 
been a number of recent mea- 
surements of the peculiar veloc- 



ity field from extant supernovae 
data (e.g., Riess, Press, & Kirsh- 
ner 1995, Zehavi et al. 1998). 

Although the TF and fundamental 
plane relations are less accurate than 
surface brightness fluctuations or Su- 
pernovae la's, they have been the 
distance indicators most widely used 
for mapping the nearby velocity field, 
mostly because they are applicable 
to the largest numbers of galaxies, 
and therefore can be used to de- 
scribe the velocity field in detail. A 
large number of surveys of nearby 
galaxies have used these techniques 
(SW; Postman 1995; Strauss 1997; 
Willick 1998). Willick et al. (1995, 
1996, 1997a) have combined the data 
from a number of these, and put 
them on a common footing, to create 
the Mark III catalog of peculiar ve- 
locities, which covers much of the sky 
to a depth of roughly 6000 km s — 1 . 
Haynes et al. (1998) have carried out 
a TF survey of spiral galaxies with 
somewhat poorer sampling, but with 
fewer difficulties concerning match- 
ing disparate data sets, to create 
their SFI data set. We discuss the 
results of these two surveys in some 
detail below. Figure [l] shows the pe- 
culiar velocity field of the Mark III 
dataset, as projected onto the Super- 
galactic Plane. Dots are drawn at the 
measured distance of the galaxy, and 
the line ends at its measured redshift 
(in the CMB frame); thus the length 
of the line is the measured radial pe- 
culiar velocity. To reduce errors, and 
the clutter in the diagram, neigh- 
boring galaxies have been grouped 
following an algorithm described in 
Willick et al. (1996). 
There is clearly 

a tremendous amount of information 
in a plot like this, but the data shown 
are noisy and sparse, and careful 
thought is needed to use it in a sta- 
tistically rigorous way. One statistic 
that one would like to measure is the 
coherence scale of the velocity field. 
That is, one can measure the mean 
bulk flow within a sphere centered on 
the Local Group; the Cosmological 
Principle implies that within large 
enough a sphere, this flow should be 
zero. The question, of course, is on 
what scale this is the case. We will 
not review this rather controversial 
subject here, other than to point out 
that the situation right now is par- 
ticularly confused (and therefore ex- 
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Mark III Velocity Field in Supergalactic Plane; CMB Frame 
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Figure 1. The Mark III peculiar velocities of all galaxies within 22.5° of 
the Supergalactic plane. The point is drawn at the measured distance of the 
galaxy, while the line is drawn to its redshift, in the CMB rest frame. Pos- 
itive peculiar velocities are drawn with solid points and solid lines, while 
negative peculiar velocities use open points and dashed lines. Points repre- 
senting groups of more than three galaxies are drawn somewhat larger. Some 
of the large superclusters in the nearby universe are marked. 



citing!): the Lauer-Postman (1994) 
measurement of bulk flow of order 
700 km s" 1 on scales of 15,000 
km s — 1 remains unconfirmed, al- 
though despite much scrutiny, no- 
body has found any fault at all 
with the data or the analysis. Mean- 
while, the SMAC survey is find- 
ing a bulk flow of a similar ampli- 
tude on similar scales, although in a 
quite different direction from Lauer- 
Postman (Smith and Hudson, these 
proceedings). On the other hand, 
Giovanelli et al. (1998a, b) are finding 
impressive evidence from their field 
and cluster TF survey for negligible 
bulk flows within a sphere of radius 
5000 km s — 1 centered on the Lo- 
cal Group. As Strauss (1997) empha- 
sizes, a clean answer to this question 



will have much to say about the mass 
power spectrum on large scales. 
In this review, we would like to focus 
on a different subject, the compari- 
son of the peculiar velocity data with 
the galaxy density field via equa- 
tion (H), and its implications for bias 
and Q. 



3 The Relative Distribution of 
Galaxies and Dark Matter 

Equation fc), or its differential form, 
equation (p), make reference to the 
mass density fluctuation field <5. We 
would like to compare the observed 
peculiar velocity field with the ob- 
served galaxy density field, as deter- 
mined, e.g., from redshift surveys of 
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IRAS galaxies (Fisher et al. 1995; 
Yahil et al. 1991; Branching these 
proceedings). This requires a model 
for the relation between the galaxy 
and dark matter distributions. 
The density fields p of galaxies and 
dark matter of course differ by their 
relative contributions to the overall 
density of the universe, so we work 
with the density fluctuation field 8 
(equation^), smoothed on some scale 
R. 

The simplest model, which was im- 
plicitly assumed in the early days of 
the subject of large-scale structure, 
was that the two density fluctuation 
fields were identical: 

'-'galaxies '-'dark matter* ('■O 

In 1984, Kaiser wrote a tremendously 
important paper, pointing out that 
if clusters of galaxies form preferen- 
tially in denser regions of space, their 
clustering would be stronger, or bi- 
ased, relative to the galaxies. It was 
quickly realized that the same state- 
ment could hold for galaxies relative 
to dark matter. Although Kaiser's 
original formulation really only made 
statements about the strength of the 
two-point correlation function, it was 
quickly made more specific with the 
so-called linear bias model: 

'-'galaxies ^ ^dark matter* (8) 

This became a quite important idea, 
as it gave modelers a crucial extra 
degree of freedom with which to fit 
cosmological theories to the increas- 
ingly stringent constraints from the- 
ory. For example, Davis et al. (1985) 
showed that the standard Q = 1 Cold 
Dark Matter model could not be rec- 
onciled with existing data without 
invoking an appreciable bias, b f-x 2.5. 
Of course, this value of b was found 
to be in strong disagreement with 
that implied by the COBE normal- 
ization, which is one of many ways 
of expressing the non- viability of this 
model. 

In addition to giving the theorist 
an extra parameter to play with, 
Kaiser's bias idea made us realize 
that there was no reason a pri- 
ori to expect that the distribution 
of galaxies and dark matter trace 
one another perfectly. The process 
of galaxy formation is complicated, 



as the various detailed models pre- 
sented at this conference by Kauff- 
mann, Lacey, and others indicate. 
Therefore, the linear proportionality 
in equation (^), with a universal bias 
parameter b that is independent of 
position, scale, and time, is almost 
certainly an oversimplification. We 
could more generally write down a 
generic relation between the two den- 
sity fluctuation fields: 

"galaxies = ■T^dark matter) 

~r e dark matter j (9) 

where J 7 is a general, nonlinear func- 
tion, which may depend on smooth- 
ing scale R, and e is a "ran- 
dom" variable, e contains all the 
physics of galaxy formation which 
does not depend on the local density 
<5dark matter* It turns out that even 
with the freedom introduced by this 
general form, one can still say use- 
ful things about the dark matter dis- 
tribution (Dekel & Lahav 1998), and 
we will come back at the end of this 
review to a discussion of what the 
properties of T and e might be. But 
for the moment, let us stick with the 
simplifying assumption of linear bi- 
asing, equation (^). In this case, com- 
parison of the peculiar velocity field 
with the galaxy density fluctuation 
field via equation (^) yields: 

qO.6 

V ■ v(r) = — <5 ga i a xics(r), (10) 

and similarly for equation (|j|). Thus 
direct comparison of peculiar veloc- 
ity and galaxy density fields allows 
us to constrain not f2 directly, but 
rather the degenerate combination 
p = Q a - 6 /b. 

4 Comparing Peculiar Veloci- 
ties with the Galaxy Density 
Field 

As indicated above, the noisiness of 
existing peculiar velocity data makes 
analyses using equation (|l(]), or its 
integral counterpart, quite difficult, 
requiring much care to control sub- 
tle systematic effects. The whole field 
is reviewed in SW; we here describe 
some of the recent approaches to the 
problem, emphasizing where their re- 
sults diverge. 
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4.1 Density-Density 
Comparison 

We would like to use equation ( |ic| ) 
to make the comparison between the 
velocity and density fields. Unfor- 
tunately, our peculiar velocity data 
give only a noisy realization of just 
the radial component of the velocity 
field at a sparse and inhomogeneous 
set of points. Dekel et al. (1990) and 
Bertschinger et al. (1990) have de- 
veloped a method, called POTENT 
(subsequently refined by Dekel 1994; 
Dekel et al. 1998) to get around this 
problem. The fundamental insight 
comes from the assumption (well- 
justified on large scales) that the ve- 
locity field exhibits potential flow, so 
that it may be expressed as the gra- 
dient of a potential <E>: 

v(r) = -V*(r). (11) 

In this case, the radial component 
of the velocity field alone contains 
enough information to determine <!>. 
In particular, integrating the ob- 
served radial component of the ve- 
locity field along radial rays yields 
<3>, and the gradient of the resulting 
field yields the full three-dimensional 
velocity field. A further divergence 
yields the Quantit y o n the left-hand- 
side of equation (DXl), which can be 
compared directly with the observed 
^galaxies from a redshift survey. 
This program has been carried out 
recently by two groups. Hudson et 
al. (1995) used an early version of the 
quantity V ■ v derived from the Mark 
III data to compare with Hudson's 
(1993) reconstruction of the density 
field of optically selected galaxies. 
They found good qualitative agree- 
ment between the two fields, and de- 
rived a value /3 op ti ca ; = 0.74 ± 0.13, 
where the optical subscript refers to 
the bias value for optically selected 
galaxies. However, they could not 
fully explain the observed scatter in 
their V ■ v — 8 comparison, given the 
known errors in each of those quan- 
tities separately, and were forced to 
invoke additional sources of scatter. 
More recently, Sigad et al. (1998) 
used the latest version of the PO- 
TENT density field to compare with 
the density field of IRAS galaxies. 
The comparison of the two density 
fields in the Supergalactic Plane is 
shown in Figure Gl The effective 



smoothing of the two fields is a 1200 
km s — 1 Gaussian. As the peculiar ve- 
locity data get sparser and noisier, 
the systematic errors in V-v get more 
difficult to control; the heavy solid 
line that outlines a shape roughly 
like the subcontinent of India repre- 
sents that region within which these 
errors are well-understood and fi- 
nite, as determined from extensive 
Monte-Carlo experiments. The qual- 
itative agreement between the two 
maps is impressive. We do a lin- 
ear scaling of one relative to the 
other (see below), take the differ- 
ence, and divide by the estimated 
errors; the result is the lower right- 
hand panel. The most significant dis- 
crepancy is found in the Zone of 
Avoidance (corresponding approxi- 
mately to SGY = in this pro- 
jection), near the Great Attractor, 
at SGX ft -40/i- 1 Mpc,S'Gy ft 
0. The peculiar velocity data are of 
course very sparse there, and the 
POTENT errors are therefore corre- 
spondingly high (causing the pinch- 
ing in the waist of the region within 
which the POTENT data are to 
be most trusted). Nevertheless, the 
discrepancy here remains a serious 
area of concern, and requires further 
work. 

If we restrict ourselves to the region 
of space within which we claim to un- 
derstand our errors (and in particu- 
lar, stay away from the Kashmir^] of 
the discrepant region near the Great 
Attractor), and do a linear regression 
of V ■ on Sjras on a grid with 
500 km s — 1 spacing (taking into ac- 
count the errors in both dimensions), 
we find a slope of 0.89. More impor- 
tantly, we can define a statistic S like 
a x 2 /dof, which relates the residuals 
from the fit to the estimated errors 
per point; we find a value S = 1.06. 
We have carried out Monte-Carlo ex- 
periments which include all relevant 
selection effects and sources of error 
to check the expectation values of 
these quantities. Our observed value 
of S falls right in the range of those 
found by Monte-Carlo experiments, 
thereby allowing us to conclude that 

3 Thanks to Luiz da Costa for this 
terminology! 

4 Actually, the analysis uses a non- 
linear extension of this linear-theory 
expression; see Sigad et al. (1998) for 
details. 
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Figure 2. The density fluctuation field in the Supergalactic Plane at 1200 
km s _1 smoothing, as determined from the POTENT analysis (upper left- 
hand corner) and the IRAS 1.2 Jy redshift survey (upper right). Contours are 
separated by 0.2 in 8. The lower left panel shows the difference between the 
two fields (after linear scaling). The lower right shows the difference, divided 
by the estimated errors; contours are separated by unity. The heavy contour 
marks the region within which the POTENT errors are well-understood. 



• we understand our errors, and 

• the data are consistent with the 
null hypotheses of gravitational 
instability theory and linear bias- 
ing. 

The Monte-Carlo experiments were 
based on an JV-body simulation with 
an effective f3 = 1; we found that the 
method delivers an unbiased result, 
with an error of order 0.12, so we con- 
clude further that 

Piras = 0.89 ±0.12. (12) 

We are in the process of carrying 
out similar experiments with both 
Q = 1 and Q = 0.3 simulations with 
(3 = 0.3, and we find that we get un- 
biased results for /3 from our analy- 



sis in this case as well, so we can say 
with some confidence that the data 
are inconsistent with (3 = 0.3. 

4.2 Velocity- Velocity Compar- 
ison 

Another approach to the determina- 
tion of (3 involves the integral form 
of the velocity-density relation, equa- 
tion (j^). In particular, this equation 
allows one to make a model for the 
expected velocity field, given the ob- 
served galaxy density field from a 
redshift survey; this may then be 
compared to observed peculiar ve- 
locities. In practice, one can avoid 
the various nasty Malmquist biases 
that would otherwise come up in 
this game by phrasing the prob- 
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lem slightly differently: establishing 
the TF relation requires knowledge 
of the distances of galaxies to de- 
termine absolute magnitudes. Given 
measured redshifts, this requires a 
model for the velocity field. The 
correct model for the velocity field 
is that which minimizes scatter in 
the TF relation (or, in somewhat 
more sophisticated approaches, max- 
imizes the likelihood of the TF ob- 
servables). One important advantage 
of this approach is that it works 
directly with the observables them- 
selves, and therefore avoids the great 
deal of data massaging that the PO- 
TENT analysis requires. 
This basic approach has been car- 
ried out by a number of workers. 
Shaya, Peebles, & Tully (1995) used 
an action principle extension of lin- 
ear theory to generate a velocity 
field model from Tully's (1987) cat- 
alog of nearby galaxies, weighted by 
the blue light of each galaxy. Min- 
imizing the scatter in the inverse 
form of the TF relation for a some- 
what massaged form of the Aaron- 
son et al. (1982) data, yielded a value 
floptical = 0.34 ± 0.130. Their ef- 
fective smoothing length was quite 
small, only a few Mpc. 
Davis, Nusser, & Willick (1996) did 
the clever trick of expressing both 
the observed (from Mark III) and 
modeled (from IRAS) velocity fields 
in terms of the same spherical- 
harmonic-based orthonormal expan- 
sion. They then developed a formal- 
ism that allowed them to express 
the inverse TF relation in terms of 
this expansion. This has the great 
advantage of smoothing the data 
on large scales to ensure linearity, 
and guarantees that the two fields 
are smoothed on equivalent scales. 
To their surprise, they found that 
within a distance of 5000 or 6000 
km s _1 , the two velocity fields dis- 
agreed in even the lower-order mul- 
tipoles. They suggested that there 
are possible systematic errors in the 
matching of the different peculiar ve- 
locity datasets which make up the 
Mark III catalog, and therefore they 
do not quote a value of ft. 
da Costa et al. (1998) used the 
same analysis as Davis et al. (1996), 

5 They expressed their results in 
terms of Q, and explicitly assumed 
b = 1 for their sample. 



now using the SFI TF sample of 
Haynes et al. (1998), which should 
be less subject to the possible match- 
ing problems of the Mark III. They 
found very good qualitative agree- 
ment between the velocity fields, and 
concluded that ftmAS = 0-6 ± 0.1. 
The final TF scatter was unfortu- 
nately appreciably larger than ex- 
pected a priori, perhaps due to com- 
ponents of the velocity field on scales 
smaller than they were modeling. 
Finally, Willick et al. (1997b) and 
Willick & Strauss (1998) carried 
out a rigorous likelihood analysis 
of the Mark III TF data (termed 
VELMOD). Unlike the other ap- 
proaches, it allows one to explicitly 
take into account the presence of re- 
gions in which the redshift-distance 
relation along a given line of sight 
is non-monotonic (triple-valued re- 
gions), the scatter due to small-scale 
unmodeled components of the veloc- 
ity field, variations of the TF scatter 
with luminosity, variations in the TF 
calibration from one dataset to an- 
other, and other systematic effects. 
VELMOD requires a full modeling of 
the velocity field, and thus the effec- 
tive smoothing of the analysis was 
pushed to as small scales as possi- 
ble; 400 km s — 1 was used in prac- 
tice. The Mark III and IRAS data 
were used here; they were found to 
be consistent with the linear biasing 
and gravitational instability model. 
The formal analysis gives a value 
ftlRAS = 0-50 ± 0.04, where the im- 
pressively small statistical error bar 
is confirmed by Monte-Carlo experi- 
ments. Interestingly, this value is in 
excellent agreement with the value 
above for optical galaxies by Shaya 
et al. (1995), when one takes into ac- 
count the relative bias of optical and 
IRAS galaxies of 1.4 (e.g., Hermit 
et al. 1996). The VELMOD analy- 
sis confirms the suspicion of Davis et 
al. (1996) of systematic errors in the 
matching of the datasets that make 
up Mark III; the best-fit TF rela- 
tions differ systematically from those 
assumed in Davis et al. (1996), in 
such a way to explain the discrepan- 
cies the latter found in their analy- 
sis. With the freedom to fit the TF 
relation, Willick & Strauss (1998) 
found excellent agreement between 
the Mark III and IRAS datasets. 
So the velocity-velocity comparisons 
carried out here are in substan- 
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tial agreement with one another. 
The density-density comparison us- 
ing the POTENT map gives a value 
of Pi has i using the same data, that 
differs by several sigma. Those who 
would argue for Q = 1 will be buoyed 
by one set of analyses, while those 
who would argue for Q = 0.3 will be 
cheered by the other. Why the dis- 
crepancy? Several possibilities come 
to mind: 

• The systematic sampling errors 
in the Mark III propagate into the 
POTENT analysis to bias the de- 
termination of Pi has- However, 
preliminary tests show that the 
POTENT density field is actually 
quite robust to these problems. 
The POTENT density field is a 
local quantity; it should be insen- 
sitive to these global problems in 
the velocity field. 

• Both the POTENT and 
VELMOD codes have been ex- 
tensively tested with Monte-Carlo 
experiments based on PM simula- 
tions of the local universe, follow- 
ing Kolatt et al. (1996). It is possi- 
ble that the finite force resolution 
of the code is underestimating 
non-linear effects on small scales, 
which may have unknown effects 
especially on the VELMOD anal- 
ysis, with its appreciably smaller 
smoothing scale. We are planning 
simulations at higher resolution 
to test this possibility. Willick & 
Strauss (1998) show that their re- 
sults are in fact quite robust to a 
substantial increase in smoothing 
scale. 

• Again, VELMOD and PO- 
TENT use quite different smooth- 
ing scales. If the effective bias is 
a function of scale, one expects 
the value of P in the two anal- 
yses to differ. Moreover, stochas- 
ticity in the galaxy-mass relation 
(i.e., the e term in equation |jj) 
can cause systematic errors in the 
determination of P (Dekel & La- 
hav 1998). Considering only sec- 
ond moments of the galaxy-mass 
relation, the regression of V ■ v on 
<5 g alaxies nas a slope of rD, ' 6 /b, 
where b = cr g /cr is the ratio of 
the variances in the galaxy and 
mass density fields and and r = 
(8g8)/o- g a is the correlation coef- 
ficient. Note that VELMOD, be- 
ing a velocity-velocity compari- 



son, uses an integral of the density 
field over a large range of scales, 
so it is more difficult to under- 
stand the effect of stochasticity 
and scale-dependent bias on these 
results. We plan to use realistic 
bias models (see below) to cre- 
ate simulations of their effect on 
VELMOD. 

4.3 Other Approaches to P and 
f! from Peculiar Velocities 

There are a number of other ways 
to constrain P and Q from pecu- 
liar velocity data and equation (JsJ) - 
The redshifts of galaxies differ from 
the distances by the radial compo- 
nent of the peculiar velocity, equa- 
tion Because only the radial 
component of the inferred three- 
dimensional position of a galaxy in 
redshift space is affected, the cluster- 
ing of galaxies acquires a systematic 
radial anisotropy. Consider a cluster 
of galaxies. On small scales, the virial 
motions within the cluster cause 
what is a compact structure in real 
space to appear stretched out along 
the line of sight in redshift space, 
thus reducing the apparent strength 
of the clustering. On large scales, 
coherent infall of galaxies towards 
this overdensity makes galaxies ap- 
pear closer to the cluster in redshift 
space than in real space, enhanc- 
ing the clustering strength. Hamilton 
(1998) has written a comprehensive 
review of attempts to measure this 
large-scale effect from redshift sur- 
veys; its amplitude can be related to 
the strength of clustering using lin- 
ear perturbation theory, yielding an 
estimate of p. His summary of pub- 
lished analyses is Piras = 0.77 ± 
0.22, but individual measurements 
are noisy, with large errors. More 
importantly, measurements of the 
anisotropy of the correlation func- 
tion or power spectrum show sys- 
tematically lower values of P than 
do analyses measuring rcdshift-space 
distortions in spherical harmonics, 
reminiscent of the density-density 
vs. velocity-velocity dichotomy in re- 
sults described above. 
Dekel & Rees (1994) realized that 
one can put a firm lower limit on 
f2, independent of P, from peculiar 
velocity data alone. By definition, 
all regions of space have 8 > — 1. 
In linear theory, we can write —V ■ 
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v = n os 8 > -O a6 . Therefore, the 
lowest observed value of — V • v al- 
lows us to put a lower limit on fi. 
Using the Sculptor Void (SGX = 
-20/i- 1 Mpc,SGy = -40/i _1 Mpc) 
in the Mark III POTENT maps 
yields a limit Q > 0.3 at the 2.4 a 
confidence level. This result needs 
further checking with simulations a 
la Kolatt et al. (1996), and with the 
latest versions of the Mark III and 
SFI data. 

One can use the POTENT maps in 
another way to put constraints on 
O. In inflationary models, one ex- 
pects the initial one-point density 
distribution function to be Gaus- 
sian, while non-linear gravitational 
growth causes systematic deviations 
from Gaussianity. Nusser & Dekel 
(1993) have developed a "time ma- 
chine" with which they can correct 
the quantity V • v for non-linearities, 
to regenerate the original linear den- 
sity field, and therefore its distri- 
bution function. This time machine 
depends on O; they find that with 
= 1, the initial distribution func- 
tion is beautifully Gaussian, while it 
is far from Gaussian for O = 0.3. 
They put a limit Q > 0.3 at Acr con- 
fidence with an early version of the 
Mark III data, under the assump- 
tion of Gaussian initial conditions. 
Again, this needs checking with mod- 
ern data, with special attention paid 
to systematics, which may effect the 
tails of the distribution. 
Finally, one can try to measure the 
power spectrum of density fluctua- 
tions that gives rise to the peculiar 
velocity field that one sees. Express- 
ing equation (g) in Fourier space, 

ik ■ v(k) = O ' 6 <5(k), (13) 

one sees that the velocity field is 
more sensitive to large-scale waves 
(i.e., smaller k) than is the density 
field. Therefore, the velocity field is 
a particularly effective way of mea- 
suring the power spectrum on large 
scales (cf. , the discussion in Strauss 
1997). Kolatt & Dekel (1997) have 
measured the power spectrum of 
V ■ v directly, while Zaroubi et al. 
(1997) and Freudling et al. (1998; 
see also Zehavi, these proceedings) 
have used a likelihood technique to 
fit the "raw" peculiar velocity data 
for the power spectrum. Of course, 
the quantity they end up constrain- 



ing is not P(k) directly, but rather 
1,2 P(fc). Comparing with models, 
and incorporating constraints from 
COBE and redshift surveys, then 
leads to a constraint on crgO 0,6 , or 
equivalently (3. These various papers 
are substantially in agreement; the 
latter paper yields Piras = 1.2±0.2. 

5 The Complications of Bias 

As we've hinted at above, one 
possible explanation of these dis- 
parate results is that our simplis- 
tic model of linear, deterministic, 
scale-independent bias is not valid. 
The analyses above differ substan- 
tially in their effective smoothing 
scales as well as in the morphol- 
ogy of the galaxies which they ob- 
serve. In addition, the methods they 
use depend differently on the de- 
gree of stochasticity in the galaxy- 
mass relationship (Dekel & La- 
hav 1998). While these issues are 
unimportant if linear bias holds, a 
nonlinear, stochastic, morphology- 
dependent bias could easily cause the 
results of these analyses to differ. To 
understand how a realistic bias af- 
fects matters, Blanton et al. (1998) 
have examined the relative distribu- 
tion of galaxies and dark matter in 
a simulation which models both the 
gravitational physics of dark matter 
and the gas physics of the baryons. 
The simulation handles star forma- 
tion by converting gas into collision- 
less particles in regions with infalling 
gas, with cooling times below the 
local dynamical time, and masses 
above the Jeans mass. One can then 
look at the relationship between the 
density field of these collisionless par- 
ticles (which we take as a proxy for 
the galaxy density field) and the dark 
matter density field. 
The resulting bias relationship is 
nonlinear, stochastic, and is a strong 
function of galaxy age. These prop- 
erties are revealed in Figure ^, which 
shows as a greyscale the conditional 
probability P(l + 5 g \l + 8) and as 
the solid line the conditional mean 
(1 + 8g\l + <5), where all quantities 
are defined with a top hat filter of 
radius 1 h -1 Mpc. Each panel shows 
the results at z = for galaxies 
formed at different epochs, as la- 
beled. Note that the oldest galax- 
ies are the cleanest tracers of the 
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Figure 3. Galaxy mass density as a function of dark matter density for each 
age quartile, at 1 /i -1 Mpc radius top hat smoothing. Each panel lists the 
range of formation redshifts included. Shading is a logarithmic stretch of the 
conditional probability P(l + 5 9 |1 + S). Solid lines indicate (1 + 5 9 |1 + 5); 
dotted lines indicate the la deviation from the mean. 



dark matter distribution: the scatter 
around the mean galaxy-mass rela- 
tion is small. However, the youngest 
galaxies show a nonlinear, even non- 
monotonic, relation with the dark 
matter; they are under-represented 
in the very densest regions of the 
dark matter map (reminiscent of spi- 
rals in the cores of clusters, although 
in the real universe clusters still rep- 
resent appreciable overdensities in 
the distribution of late-type galaxies; 
Strauss et al. 1992). Also, the scatter 
around the mean density relation for 
the youngest galaxies is large. 
In these simulations, the relation- 
ship between galaxies and mass also 
depends on scale. In Figure we 
show the bias 6 = c g /a calculated 
on various scales. The obvious scale- 
dependence of b is due to the de- 
pendence of the galaxy formation 
process on temperature. The tem- 
perature sets the local Jeans mass, 
which partly determines whether 
star-formation occurs: the higher the 
temperature, the greater the over- 
density needed to form stars. On 
small scales the temperature is pro- 
portional to the gravitational poten- 
tial <fi. Note that in Fourier space, 



4>(k) oc S(k)/k 2 . For high k, then, 
there is little power in the poten- 
tial or temperature fields; i.e. these 
fields are smoother than the density 
field. Thus, temperature correlates 
over large scales; furthermore, on 
these large scales it correlates with 
density as well. Thus the dependence 
of galaxy formation on temperature 
can couple the galaxy density on 
small scales with the dark matter 
density on larger scales. As Blan- 
ton et al. (1998) show, this coupling 
causes scale-dependence of the bias 
relation. The dependence of galaxy 
formation on local gas temperature is 
likely to be important in any galaxy 
formation scenario; thus, this scale- 
dependence may be generic. 
The work ahead is evaluate the 
consequences of this complicated 
bias relationship for statistical mea- 
sures of large-scale structure. The 
data definitely allow stochasticity; 
in the POTENT-JiMS comparison 
of Sigad et al. (1998), although the 
value of x 2 /dof is consistent with 
deterministic bias, it is also consis- 
tent with the inclusion of the rms 
value of e found in the simulations 
at 1200 km s _1 smoothing. More- 
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Figure 4. The bias b(R) = a g (R)/cr(R), where R refers to the top hat 
smoothing radius. Solid line indicates all galaxies. Dashed lines indicate each 
age quartile, with range of formation redshifts listed. Note the strong scale- 
dependence, and that old galaxies are more biased than young. 



over, nonlinear relative bias between 
galaxies of different types is unam- 
biguously observed, in the form of 
the morphology-density relation (cf. 
the discussion in Strauss & Blanton 
1998). We are in the process of car- 
rying out simulations using realistic 
bias laws, to determine how it may 
skew the results of the various anal- 
yses described above. 
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